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ABSTRACT 

We present HST/NICMOS Paschen a (Paa) images and low and high resolution IRS spectra of 
photoevaporating disk-tail systems originally detected at 24/im near O stars. We find no Paa emis- 
sion in any of the systems. The resulting upper limits correspond to about 2—3 x 10~^ Mq of mass 
in hydrogen in the tails suggesting that the gas is severely depleted. The IRAC data and the low reso- 
lution 5—12 /xm IRS spectra provide evidence for an inner disk while high resolution long wavelength 
(14—30 /xm) IRS spectra confirm the presence of a gas free "tail" that consists of ^ 0.01 to ~ 1 /im 
dust grains originating in the outer parts of the circumstellar disks. Overall our observations support 
theoretical predictions in which photoevaporation removes the gas relatively quickly (< 10^ yrs) from 
the outer region of a protoplanetary disk but leaves an inner more robust and possibly gas-rich disk 
component of radius 5-10 AU. With the gas gone, larger solid bodies in the outer disk can experience 
a high rate of collisions and produce elevated amounts of dust. This dust is being stripped from the 
system by the photon pressure of the O star to form a gas-free dusty tail. 
Subject headings: circumstellar matter - planetary systems: protoplanetary disks - stars: formation 
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1. INTRODUCTION 

The timescale at which a protoplanetary disk loses its 
primordial gas, which dominates the total disk mass, is a 
crucial factor regarding the outcome of the planet form- 
ing process in the vicinity of high mass stars. The pho- 
toevaporation of protoplanetary disks by external UV 
radiation in an environment dominated by high mass O- 
stars is a particularly potent mechanism for removing 
disk material. The timescale of photoevaporation is gov- 
erned by the amount of UV radiation from the nearby O- 
stars, which significantly erodes the protoplanetary disk 
in less than a million years. Studies of disk erosion agree 
that the outer part of the disk evaporates qu i ckly; how- 
ever according to, e .g., lHollenbach fc Adami (J2004D and 
lAdams et al.l ()2004[ ) the erosion of the disk slows down 
signi ficantly when the size reaches 10 AU. Other au thors 
(e.g. lJohnstone et al.|[l998l : iMatsuvama et al.ll2003D have 
found that the photoevaporation process can shrink the 
radius of the disk below 1 AU within 10^ years and it 
can even destroy the whole disk when it is coupled with 
internal viscous evolution. Therefore, the presence of 
high mass stars might be lethal for forming planets since 
the EUV and FUV radiation of a nearby 0-star can dis- 
rupt the protoplanetary disk on a short time scale (a 
couple of times 10^ yr) and stop the planet formation 
process. On the other hand the radiation might trigger 
the formation of rocky planets by accelerating the dust 
settling and grain growth in the midplane of the disk 
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(iThrooD fc Ballvl[2005h . 

iBalog et al.l (|2006D reported the discovery of three 
comet like structures near different 0-type stars (with 
the comet-like tails pointing away from the central high 
mass star) in 24 iivn Spzteer/MIPS images. A dopting a 
radia tion-pressure driven outflow disk model (jSu et al.l 
l2005l they concluded that these objects are photoevapo- 
rating protoplanetary disks where dust entrained in the 
photoevaporative flow is responsible for the mid-infrared 
tails detected at 24 /xm. The tails extend up to 0.2 pc 
from the evaporation working surface, so they are conve- 
nient for detailed study. To investigate this phenomenon 
further and to get a better understanding of the nature 
of these objects, we have searched these three systems 
for the gas expected to be a signiflcant component of the 
flow. We used Paa as an ideal tracer of the photoevap- 
oration process, since it is the strongest emission line in 
the near infrared and is not strongly affected by extinc- 
tion. We also report mid-IR spectra to provide a better 
understanding of the physical processes in the disk and 
to help determine the properties of the dust. 

2. OBSERVATIONS 

2.1. HST/NICMOS 

We used HST/NICMOS to search for ionized hydro- 
gen, predicted to be an abundant p roduct of photoe- 
vapo ration in a protoplanetary disk (jRichling fc Yorka 
I2000D . We chose the fllters F187N and F190N (Paa line 
and a nearby lineless continuum) because we expected 
to detect the strongest emission in these filters. Our 
targets wer e the three objects detected at 24 ^m by 
iBalog et al.l ()2006[) . All these objects have a cometary 
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Fig. 1. — The images of the photoevaporating source in 
NGC 2244. Upper left: 24 /im; upper right Pact; lower right: Paa 
continuum; lower left continuum subtracted Paa. The covered area 
is 19" X 19" (0.138 X 0.138 pc in the distance of NGC2244) in all 
four panels. 



structure in the 24 /im Spitzer/MIPS images with the 
"comets' tail" pointing way from a nearby O star (see 
Fig. 1 of lBalog et aH ^QM))- 

The exposure time for both the hne and the contin- 
uum filters was 272 s, 352 s and 671 s in the case of ob- 
jects in NGC 2244, NGC 2264 and IG 1396 respectively. 
We employed scale factors of 3.2641699x10^5 Jy s DN^^ 
for F190N and 3.2329899x10^5 Jy s DN"i for F187N 
in the case of NGC 2264, 3.2536798x10-5 Jy s DN^i 
for F190N and 3.1735399x10-^ jy s DN^i for F187N in 
the case of NGC 2244 and 3.2641699x10^5 Jy s DN^i 
for F190N and 

3.2329899x10-5 Jy s DN"! for F187N in the case of 
IG 1396 to convert the original frames to calibrated flux 



images. 

As an example Fig. [1] shows the source in NGC 2244 
at 24 iim, at 1.87 /im ( Paa line) and in the continuum 
subtracted Paa images. The covered area is 19" x 19" 
(0.138 x 0.138 pc in the distance of NGC2244). The 
other two sources behave similarly. It is clear that the 
tail structure is not visible in the Paa images even after 
continuum subtraction. However in all cases there is an 
unresolved point source in the comets' head, which con- 
firms that the photoevaporating source is indeed a star 
and not an extended globule. 

We performed aperture photometry on the continuum 
images following the prescription in the NICMOS data 
handbook to estimate the continuum flux of the central 
object. In two cases (NGC 2264 and IG 1396) the NIC- 
MOS fluxes are in good agreement with fluxes calculated 
from the 2MASS data (See Figs. [3and[SD. For NGC 2244 
the nearby O star strongly affects the 2MASS magni- 
tudes so the fluxes calculated from the 2MASS data are 
much higher than the flux we obtained using NICMOS. 
However VLT/VIMOS R and I fluxes (Robert King pri- 
vate communication) confirm that our NICMOS mea- 
surements are accurate and the 2MASS data are about 
1-2 magnitudes higher than the real values. Therefore 
we ignore the 2MASS data for this source and use only 
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Fig. 2. — Spectrum of the head of the photoevaporating source 
in IC 1396 (light blue line: low resolution; magenta line: smoothed 
high resolution). The green line represents the continuum of a star 
with stellar temperature of 4000 K and blue line represents the 
SED of our dust tail model while the red line is the sum of the star 
and tail models. Black dots show the fluxes of the central point 
source in different photometric bands, orange dots show the fluxes 
of the tail (we give 3cr upper limits where the tail is not detected) 

the R, I and NICMOS 1.9 ^m data to estimate the SED 
of the central object (Fig. [4|). The reason why this prob- 
lem does not affect the other two objects is that they are 
about five times farther (in angular distance) from the O 
star than the source in NGG2244. 

2.2. Spitzer/IRS 

We used Spitzer/IRS to take low and high resolution 
mid-IR spectra of the photoevaporating sources to ex- 
plore the physical processes in the evaporating gas. All 
spectra were taken centered at the "head" positions. 
The exposure times for the four IRS modules (short-low 
1st order (SLl), short-low 2nd order (SL2), short-high 
(SH), and long- high (LH)) are given in Table 1. We ex- 
tracted the low resolution spectra from the Spitzer Sci- 
ence Center (SSC) calibrated post-bed background sub- 
tracted frames using SPIGE^. We averaged the two nod 
positions to increase the S/N. For the high resolution 
modes we subtracted a spectrum obtained at a nearby 
location before we extracted the science spectra, to elimi- 
nate contamination from the surrounding HII region. We 
also discarded the beginning and the end of each order. 
We calibrated the extracted spectra using the 5.8—24 fim 
photometry. No emission lines are apparent in the high 
resolution spectra, indicating the absence of gas. How- 
ever, the Paa observations place more stringent limits on 
the amount of gas. We smoothed the IRS high resolution 
spectra to reduce the noise, and then merged them to get 
full coverage between 10 and 30 /im. 

We show our spectra in Figs. [21 [Sj and [4] along with 
the Kurucz model spectra for the stellar photospheres 
and our model SEDs of unbound small particles. The IR 
excess is clearly visible even in the unsmoothed low S/N 
low resolution spectra. The high resolution spectra show 
a bump around 11.3 /Ltm that is not detected in the low 
resolution spectra. This feature originates from aromatic 
emission and is very common in HII regions. The reason 
we detect it at high resolution and do not detect it at 
low resolution can be explained by the difference in the 
background subtraction in the two observing modes. In 

2 Spitzer IRS Custom Extraction (SPICE) is a JAVA-based tool 
for interactively extracting Spitzer/TRS spectra. 
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TABLE 1 
The Spitzer/IRS observations of the photoevaporating tails. 



Cluster 


J2000 

[h:m:s] [d:m:s] 


SLl 

[s] 


Integration 
SL2 

[s] 


SH 

[s] 


LH 


NGC 2244 
IC 1396 
NGC 2264 


06:31:54.68 04:56:25.0 
21:38:57.09 57:30.46.5 
06:41:01.92 09:52:39.0 


18.1 
146.8 
146.8 


73.4 
146.8 
146.8 


25.1 
487.5 
125.8 


31.5 
2176.4 
457.6 
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Fig. 3.— Same as in Fig[2] for the source in NGC 2264 
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Fig. 4. — Same as in Fig[2] for the source in NGC 2244. We note 
that the low S/N in the short wavelength portion of the high res- 
olution spectrum resulted in two spurious absorption-like features 
in the final smoothed spectrum. These features are not real so we 
do not discuss them any further. 

the short-low (SL) mode we subtract the two nod posi- 
tions from each other to ehniinate the background while 
in high resolution mode a different spectrum taken at a 
nearby location was used for the same purpose. We see 
the 11.3 /im line in the latter case as a result of imper- 
fect background subtraction due to spatial variation of 
the line flux. To support this explanation we show the 
relevant part of the low resolution spectrum of the source 
in NGC 2264 (Fig. [Sjleft panel) and the same part of a 
low resolution spectrum taken at the nearby background 
location where the high resolution background spectra 
were taken (Fig. [5] right panel). The two sharp verti- 
cal lines are the spectra of our target and a nearby star 
that happened to fall on the slit while the smeared out 
horizontal line in the middle of the image is the 11.3 /xm 
aromatic feature. The facts that the aromatic feature 
is visible through the entire slit width and that it is in 
both spectra suggest that it comes from the background 




Fig. 5. — Part of the low resolution spectrum of the source 
in NGC 2264 (left panel) and the same part of a low resolution 
spectrum taken at the nearby background location where the high 
resolution background spectra were taken (right panel). The ex- 
posure times for the two spectra are exactly the same as well as 
the stretch of the images as indicated by the greyscale code at the 
bottom of each image. The direction of the dispersion is along 
the vertical axis while the horizontal axis represents the spatial 
direction along the slit. 

rather than from the source itself. It is also clear that 
the feature is much weaker at the background location 
(right panel) than at the source position (left panel). 

2.3. Archival Spitzer/IRAC-MIPS data 

We downloaded all the available 3.6-8 /im data from 
the Spitzer archive to see if new observations have be- 
come available since the appearance of the paper report - 
ing the discovery of these objects (jBalog et al.l 120061 ). 
Only in the case of NGC 2244 did we find new deeper 
IRAC observations (Bouwman et al., Prog ID: 30726). 
For the rest of the sources we used the data from the 
Prog ID: 37 (PI: Fazio) and Prog ID: 58 (PI: Rieke) GTO 
programs. 

We re-mosaiced all the IRAC observations using cus- 
tom IDL scripts and extracted the flux at the position 
of the "head" of each source using 2 pixel apertures. We 
estimated 3ct upper limits for the flux of the "tail" in 
the case of NGC 2264 and IC 1396. The situation for 
NGC 2244 was more complicated since the nearby O star 
affects the photometry at shorter wavelengths. As a re- 
sult, although we detect the source at 3.6 and 4.5 /im 
we cannot extract its flux. At 5.8 and 8.0 /im both the 
"head" and the "tail" are measurable in the new mosaics. 

For the 24 /im MIPS photometry we used the mo- 
saics originally generated with the DAT ([Gordon et al.l 
20051) and pubhshe d by JSicilia- A gmlar et al I (|2006f) . 
Teixeira et all ()2006 ^. and Balog et al.l (|2007l ). First we 
measured the flux of the "head" via PSF fitting using 
IDP3^ then we estimated the flux in the tail. To do so, 
we selected a rectangular region along the tail in the PSF 
subtracted image and measured the signal in this region. 
Next we moved the selected region to nearby background 
positions on both sides of the tail and averaged the sig- 
nal measured at these positions. We then subtracted the 

^ Image Display Paradigm 3 (IDP3) is an IDL (Interactive Data 
Language) package for analyzing astronomical images developed 
by the NICMOS Software Group at the University of Arizona. 
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averaged background from the signal at the position of 
the tail to get the total flux in the tail. The results of 
our photometry are plotted in Figs. [H [21 and 3] (black 
dots: head, orange dots: tail). 

3. RESULTS 

3.1. Calculation of Hydrogen Mass Upper Limit 

We detected no Paa emission in any of the objects; 
we will now derive upper limits to the amount of gas 
compatible with this result. 

3.1.1. Photoevaporation from an Optically Thick (to Lyman 
continuum,) Disk 

Consider a spherical (or hemispherical) cloud of hydro- 
gen with a radius r, located at a distance d from an O 
star that has a Lyman continuum luminosity Q (in 73"^). 
If the cloud is sufficiently dense, it will be optically thick 
to the Lyman continuum radiation. I n the model de- 
velop ed for the Orion Nebula proplyds (jJohnstone et al.l 
Il998| ). soft, non- ionizing UV heats the outer layers of the 
cloud resulting in a quasi-steady-state, soft-UV driven 
wind expanding with a velocity of a few kilometers per 
second. Ionizing Lyman continuum radiation forms a 
quasi-stationary ionization front (I-front) in this wind 
where the gas becomes ionized and heated to about 10^ 
K. The resulting pressure jump drives a D-type shock 
into the expanding low-velocity neutral wind that com- 
presses and decelerates the flow before it passes through 
the I-front. At the I-front, the plasma accelerates to 
about the sound speed in the ionized gas, c//, and the 
pressure gradient results in a spherically divergent flow. 
Once ionization equilibrium has been established in the 
ionized flow, and the system is in a steady state of photo- 
ablation driven mass loss, the flux of Lyman continuum 
incident on the cloud, F = Q/Aird^ is balanced by the re- 
combination rate between the I-front and the illuminat- 
ing star. For a constant velocity, spherically divergent 
flow, ne(r) — nj{ri/rY where nj is the plasma den- 
sity just outside (on the irradiated side) of the ionization 
front, and r/ is the ionization front radius. The total re- 
combination rate in the column between the O star and 
the I-front is then given by 



Q 

47rd2 



nl{r)aBdr 



njasri 



(1) 



where the integration limits run from rj to infinity and 
as ~ 2.6 X 10~^^{cm^s^^) is the case-B recombina- 
tion coefficient for hydrogen at 10^ K. For a disk or cir- 
cumstellar environment with a mean neutral gas density 
n{H) >> (SF/aBriY^^ inside radius r/, the electron 
density at the I-front is self-regulated by the Lyman con- 
tinuum flux to be n/. A higher electron density would 
shut off the UV, decrease the mass loss, and therefore 
lower the density. A lower electron density would result 
in a higher flux at the I-front and a greater mass-loss 
rate until the density at the I-front reached the value n/ 
. If the mean neutral gas density is lower than nj, the 
I-front would fully ionize the circumstellar environment 
which would then be optically thin to Lyman continuum. 
This simple model ignores the expected increase in ve- 
locity of the photo-ablating plasma due to the density 
and pressure gradient. In a full treatment of the hydro- 
dynamics, the flow velocity will increase to about 2 to 



3 xc// as the density decreases by one to two orders of 
magnitude. However this acceleration will have little ef- 
fect on the observed surface brightness of species such as 
Paa because most of the line emission is produced by the 
dense plasma at the base of the ionized flow, since line 
intensities are proportional to the emission measure. 

The emission measure is defined in the usual manner 
to be EM = J n1{l)dl where rie is the electron density 
and dl is an infinitesimal length element along the line- 
of-sight. For a hemispherical or spherical surface giving 
rise to a constant velocity wind with an r^^ density pro- 
file, the highest EM occurs along a line-of-sight that just 
grazes the ionization front. By integrating the line-of- 
sight along a tangent to the surface, it can be shown that 
EM = {7r/2[pc])n'jri where pc = 3.086x lO^^cm converts 
the units of EM from c.g.s. to cm^^pc. From ionization 
equihbrium, n^rj w SQ/ATrasd'^, independent of the ion- 
ization front radius, so that EM « 3Q/8[pc\aBd^. The 
approximate equality is used here because we ignore the 
slight acceleration of the ionized front due to the pres- 
sure gradient. Thus, the expected surface brightness of 
any emission line is, to first order, independent of the 
cloud radius and inversely dependent on the square of 
the distance from the source of ionization. 

The flux of any recombination line of hydrogen can be 
conver ted into an est imate of the emission measure. Fol- 
lowing [SBitzefl (|1978( ). the intensity of a transition from 
level m to level n is given by 



Imn = / judi^ = 



f^^riin^r, 



nl{l)dl = 



Att j ^' ' (2) 

2.46 x IQ^'^hVr^na^nEM 

where amn is the total recombination into the upper level 
of the transition taking account of all recombinations into 
higher levels that cascade into the level of interest and 
of the branching ratio for transitions to lower levels. EM 
is the emission measure in units of cin~^pc. For full 
ionization, n^ — Up — n{H) is assumed. The intensities 
of various hydrogen recombination lines relative to II/3 
were tabulated by PcngcUy (1964). From his tables, the 
Case-B recombination coefficient for Paa is 043 = 4.10 x 
10~^''(cm'^s~^). Thus, the relationship between Paa flux 
in units of ergs~^cm~^arcsec~^ and EM is Fp^ = 2.51 x 
10~^^-BM where EM is in units of cm~^pc. 

The Paa detection limits for our three sources are given 
in Table 2 along with the exposure times. Table 3 lists 
the nearest ionizing source in each cluster along with its 
estimated spectral type and Lyman continuum luminos- 
ity, the projected separation between the head of each 
tail and this star, the expected flux of Lyman continuum 
radiation in photons per square centimeter per second, 
and the predicted flux of Paa emission under the assump- 
tion that each comet tail head is optically thick to Lyman 
continuum radiation. This estimate ignores the attenua- 
tion of the Lyman continuum radiation by dust. Finally, 
the last column in Table 3 gives the ratio of the obser- 
vational limit on the Paa emission derived from Table 2 
divided by the expected flux under the assumption that 
each star with a dust tail contains a hemispherical, high 
density core that is opaque to the Lyman continuum. 

In these estimates, we use the projected distance be- 
tween the nearest ionizing star and each tail, an assump- 
tion that will underestimate the true separation and over- 
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TABLE 2 
The HST/NICMOS Paq observations of the photoevaporating tails. 



Cluster 



J2000 

[h:m:s] [d:m:s] 



Integration 



-Fpo, limit 

— 9 —1 



NGC 2244 06:31:54.68 04:56:25.0 272 

IC 1396 21:38:57.09 57:30.46.5 671 

NGC 2264 06:41:01.92 09:52:39.0 352 



4.2 X 10-lS 
1.6 X 10-1^ 
2.9 X 10-1^ 



TABLE 3 
Physical properties of the environment of the photoevaporating tails. 



Cluster 



O-Star Spect. Type 



Q 



7S 



'11 



d 

fcml 



F(LyC) Fpa predicted Fpa deficit 



[^1 [ 

[ cm ^ s J I 



cni'^ s arcscc^ 



NGC 2244 HD4615 06Ve 

IC 1396 HD206267 06Ve 

NGC 2264 S Mon 07Ve 



2.2 X 10*3 

2.2 X 10*^ 

1.3 X 10*9 



3 X 10^^ 

1.05 X 10^** 

9.6 X 10" 



1.9 X 10^^ 
1.6 X 10^2 
1.1 X 10^2 



2.9 X 10 
2.3 X 10 
1.7 X 10 



-11 

-12 
-12 



1.5 X 10-5 
7.0 X 10-5 
1.7 X 10-* 



estimate the incident Lyman continuum flux. However, 
we also ignore the contributions of other OB stars in each 
HII region to the total Lyman continuum radiation field. 
This assumption will tend to result in an underestimate 
of the incident radiation field. The estimated incident 
radiation fields are probably uncertain to about a fac- 
tor of two. Thus, the predicted emission measures and 
Paa fluxes are order-of-magnitude estimates. However, 
as shown by Table 3, the predicted fluxes are many orders 
of magnitude higher than the observational limits. 

The observations indicate the 5piteer-detected tails are 
deficient in Paa emission by about a factor of one mil- 
lion compared to a simple model in which each dust tail 
contains a photo-ablating core of dense material. As dis- 
cussed above, the emission measure, and therefore the 
fluxes of hydrogen recombination lines are to first order 
independent of the radius of the ionization front. They 
depend on the incident Lyman continuum radiation field, 
and inversely as the distance between the object and 
source of ionization squared. Thus, a set of optically 
thick, neutral condensations ought to increase in Paa 
flux as the inverse square of the separation. If these dust 
tails were similar to the proplyds in the Orion Nebula, 
they would be orders of magnitude brighter than the ob- 
served limits. 

3.1.2. The Optically Thin (to Lyman Continuum) Case 

The absence of detectable hydrogen recombination line 
fluxes might be explained if the disks are optically thin 
to Lyman continuum radiation and the gas-to-dust ra- 
tio is signiflcantly lower than in the ISM. For an as- 
sumed radius of the circumstellar environment, rj — 
10'^ AU, the critical density above which the gas is 
opaque to the flux of Lyman continuum radiation is 

ni > {iF/aBTi'^^'^ 



2.8 X lO^FiYVg"^/^ where F12 is in 
-2 



units of 10^^ Lyman continuum photons s^^cm^"^ and ry, 
is in units of 10^ AU. This implies an emission measure, 
EM « 6 X l{f'Fi2 cm-^pc. In contrast, our Paa flux 
limits (Table 2) imply emission measures of EM « 1673 
cm-^pc for the Paa flux hmit from the NGC 2244 tail 
and EM « 637 cm^^pc for the Paa fluxes from the tail 
in IC 1396. Thus, for an incident Lyman continuum flux 
Fi2, the Paa deficit would be between 2 and 12 x 10^''. 



Table 3 gives the actual deficits for each tail given the 
projected separations from the nearest 0-star, the Ly- 
man continuum luminosity of the star (also in Table 3), 
and the flux limits of Table 2. 

If the circumstellar environment is optically thin to 
the Lyman continuum and has a constant density and 
radius r/, the upper limit for the Hydrogen mass is 

M « /xm//(11.18i?M[pc])^/^r^ , which implies masses 
of 1.7 x IQ-^Mq and 2.7 x IQ-^Mq for our two limit- 
ing cases of EM « 637 cm-^pc and 1673 cm^^pc re- 
spectively. These are very conservative upper bounds 
since observations show that circumstellar environments 
have flattened disk-like geometries. The emission mea- 
sure and derived density limits would be an order of 
magnitude lower for a disk geometry. The mass of 
a typical protoplanetary d isk is about 0.1-0.001 Mq 
() Andrews &: Williams! l2005l ). Our inferred gas mass up- 
per limits for the optically thin case are at least two 
orders of magnitude lower than the gas masses expected 
in a typical T Tauri disk, implying gas to dust ratios 10^ 
to lO'' times lower than in the ISM . 

3.2. The Inner Disk 



iBalog et"al] (|2006l ) modeled the grain properties and 
tail morphologies for the three evaporating disks dis- 
cussed in this paper. The models utilized optical con- 
stants for astronomical silicates and assumed a single 
grain size. It was found that the tail morphologies (faint- 
ness at 8 /im and the apparent lengths) placed strong 
constraints on the grain sizes. The minimum sizes were 
of order ~ 0.01 /im in radius - smaller grains emitted too 
strongly at 8 fini to be included in the tails in signifi- 
cant numbers. Grains above ~ 1 /im in radius tended 
to be too cold to fit the tail lengths. The spectra pre- 
dicted by these models for the objects in IC 1396 and 
NGC 2264 are shown in Figures [2] and [3] for comparison 
with our IRS spectra. The fit is satisfactory at wave- 
lengths longer than ~ 12 /im. At shorter wavelengths, 
there is significant emission from the "comet head" both 
from the stellar photosphere and from an excess bridg- 
ing from 2.4 /im to 12 /im and centered on the head. 
The small difference between the flux of the tail and our 
model SED at 24 /im is due to the fact that the slit of 
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LH module of IRS covers the "head" and some part of 
the "tail" so the spectra show the flux coming from some 
combination of the "head" and "tail" fluxes. We scaled 
both the spectra and the model SEDs to the flux com- 
ing through the slit to make the agreement between the 
model and the spectra at wavelengths longer than 12 fim 
more visible. For the source in NGC 2244, the previous 
model was modified to improve the fit (see Figure |4]). We 
used the new 8/im and 24 fim fiuxes to recalculate the 
dust mass and the SED for this system. The dust mass 
in small (0.03 /xm) grains is 3.5 x lO^^M©. We achieved 
a similarly good fit to the surface brightness distribution 
with the new mass estimates and the new data that are 
shown in Fig. 2 ofLBaloe ct al. (2006). 

To be conservative, IBalog et al.l (|2006ll computed dust 
masses assuming that the mass spectrum extended as 
a~^-^ (where a is the grain size) up to 1mm. The mass 
goes essentially as the square root of the maximum grain 
radius in the size spectrum. Anticipating the discussion 
in Section 4.2, the masses for our models assume a max- 
imum size of a = 1 fiia. They are entered in Table 4. 

The fits to tail spectra fall far short of accounting for 
the excess emission above the stellar photospheres in the 
3 to 10/xm range, apparent in both our spectra and IRAC 
photometry. We conclude that these excesses indicate 
the presence of class II sources in the "heads" of our 
objects. Further tentative evidence of the presence of a 
class II source in the "head" region is a broad 10 ^m 
silicate emission feature detected in the spectrum of the 
object in IC 1396. In the case of NGC 2244 the S/N of 
the low resolution spectrum is not enough to detect such 
a feature confidently. On the other hand in the case of 
NGC 2264 the SED of the unbound dust is rising very 
sharply between 10 and 15 ^m and this rise might mask 
the silicate feature at 10 /im. 

3.3. Gas in the Inner Disk 

The presence of an inne r compact disk agree s 
well with the models of iJohnstone et al.l (I1998D : 
HoUenbach & Adams' (2004"); 'A dams et all (|2004f ) and 
Throop fc Ballv (2005) which all predict a gas rich inner 
disk surviving the photoevaporation process for longer 
than 10^ years. The SEDs of the "comet" heads pre- 
sented in Figs. O m and m are very similar to the lower 
qu artile of T Tauri dis ks in the Taurus region presented 
bv lFurlaneraLl(|2006D . 

However the Paa places stringent limits on the gas 
content of even a compact disk. If a gas disk exists inside 
5 AU (which is the gravitational radius at 10 km/s, the 
sound-speed in ionized H with a cosmic abundance), its 
surface will be photo-ionized, forming a gravitationally 
bound (to the central star) HII region. It is ionized by the 
external O star, so its density is given by the Stromgren 
density. 



L{LyC) 



Anat 



0.5 



5 X lO^LOgSMgO-^cio 






(3) 



where, rji is the gravitational radius, r// = GM/cjj, 
where c// = 10 km/s is the sound speed in photo-ionized 
gas, M is the mass of the star in the disk in Solar mass 
units (1 Mq), L^g is the Lyman continuum luminosity 
of the O star in units of 10''^ photons/sec, and D is the 



distance between the O star and the disk in units of 1 
pc. 

Thus, a typical disk corona around a Solar mass star 
will have an electron density (ne2) of order 50,000 per 
cubic cm. The emission measure will be EM — ne2L — 
1.2 X lO^cm^^pc for L = lAU. This is a lower bound 
because it assumes that the central star makes no contri- 
bution to the Lyman continuum luminosity. This large 
EM will produce a strong Paa line that should be de- 
tected in the narrow band Paa images. 

The previous calculation applies only for a spatially 
resolved disk because the emission measure is relevant 
to the surface brightness assuming that the disk corona 
is resolved. In our observations the expected size of the 
inner disk (5-10 AU) subtends a much smaller angle than 
the resolution element. To apply the calculations to our 
data (in the case of NGC 2244) we consider an r^ - 5 AU 
radius disk, located at d = O.lpc from an 05 star. We 
assume that the distance from the Sun is D = 1500 pc, 
and that each Lyman continuum photon absorbed by the 
disk results in the emission of / = 0.1 Paa photons. 

In photo-ionization equilibrium, the Paa fiux from the 
disk is 

47rd^ ^'d 



F{Paa) 



f 



AtiD^ 



(4) 



This is the flux of intercepted LyC photons, multiplied by 
/, divided by AttD^, to give the fiux at Earth {L{LyC) — 
1049, rd = 5 AU, d = 0.1 pc, D = 1.5 kpc). The resulting 
fiux in this case is F(Paa)=1.27 x 10""'^^ erg cm~^ s~^. 
According to the NICMOS Exposure Time Calculator 
this fiux would provide a signal-to noise ratio around 1 
so a disk with a radius of 5 AU would not be detected. 
However the same calculation for a disk with radius of 
10 AU yields SNR>4 which would be detectable in the 
NICMOS images, implying an upper limit between 5 and 
10 AU for the radius of the inner gas rich part of the disk. 
This is in agreement with the photoevaporation models. 

4. DISCUSSION 
4.1. Mass Loss Rates and Amounts 

As we showed in Section 3.1, the Paa images and the 
mid-infrared spectra suggest that the gas-to-dust mass 
ratio in these systems is at least two orders of magni- 
tude lower th an the con ventional value of 100:1 that was 
assumed in B alog et al.l (2006), so the fiow is essentially 
gas free. 

The low gas-to-dust mass ratio requires changes in the 
assumptions made when calculating the mass loss rate 
of these obj e cts. I nstead of the value 10 km/s used in 
Balo g et al.l ()2006D assuming that the gas is dragging 
the dust, we have to use the terminal velocity assum- 
ing that the primary forces on the grains are radiation 
pressure from the O star and a drag force from parti- 
cl es in the surroundin g HII region. We use Eq. B.4 
of lAbergel et al.l (|2002f ) to estimate the terminal veloc- 
ity. Our calculations show that it is about an order of 
magnitu de higher (around 2 00-300 km/s) than the value 
used by iBalog etaTI ()2006l ) . However the new gas-to- 
dust mass ratio decreases the total amount of outfiowing 
material by about 2 orders of magnitude leaving us with 
a mass loss ra te about 10 t i mes lo wer than was originally 
calculated by IBalog et al.l (|2006[ ). Although the lower 
limit on the mass loss rate of the dust is an order of 
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magnitude higher (« 10^^^ — lO^^^M0/yr), the upper 
limit has not changed significantly because the change in 
the assumption of the maximum grain size (from 1mm 
to 1 /Ltm; see discussion in Section 4.2) and the changes 
in the assumption about the outflow velocity and in the 
gas-to-dust ratio cancel each other (leaving only a factor 
of three difference) resulting in mass loss upper limits of 
« 10-1° - lO-"M0/yr. See Table 4 for the values for 
e ach individual case. 

IStorzer fc HoUenbachl (|1999|) analyzed 10 proplyds and 
found evidence that the photoevaporative flow is dust 
and metal poor sugges ting that most s o lids a re settled 
in the disk mid-plane. iThroop fc Ballvl ()2005l ) modeled 
the evolution of a circumstellar disk in the presence of 
external UV radiation and suggested that as a result of 
photoevaporation the outer disk becomes completely gas 
depleted. It is possible that UV radiation has photoab- 
lated the disks for our three sources down to the grav- 
itational radius in ionized gas (~3-5 AU for 1 M0) so 
these cometary structures are produced by a highly pro- 
cessed disk that has lost most of its gas and started the 
dust regeneration process in the midplane of the disk. 
Radiation pressure from the central 0-star and dynam- 
ical processes in the circumstellar disk would then be 
expelling this dust from the disk. 

Assuming nearly radial orbits and relative velocities 
between the 0-star and the disk-bearing star equal to the 
cluster velocity dispersion we can estimate the time the 
disk-bearing stars have spent in the spher e of influence 
of the 0-star (< 0.5 pc iBalog et al.ll2007 t) in our three 
systems. Assuming also that the current distances of the 
disk-bearing stars to the 0-stars are the closest point of 
their orbit to the 0-star (giving the lower limit for the to- 
tal time spent in the 0.5 pc sphere) we find that they have 
spent about 7 — 10 x 10^ years within the sphere of influ- 
ence. We can estimate how much material they lose until 
they arrive at their current position based on the pho- 
toevaporation rate-distance relation of iRichling fc Yorkd 
dmi) {A'iph = 1.4 6 X lO-^Me yr-HTni^)"")- We 
adopted a — 1 from B alog et al.l ()2007[ ). We integrated 
the total mass loss during the time the disk-bearing star 
spends in the sphere of influence of the 0-star. Our cal- 
culation showed that in our three systems the total lost 
disk material is about 0.012-0.016 Mq. We present the 
results in Table 4. 

The total lost mass is sufficiently large that a typ- 
ical protoplanetary disk {M ~ 10^^ to lO-^M© 
[Andrews fc Williams! |2005[) would be totally photoevap- 
orated during the time it spends in the close vicinity of 
the 0-star. We also calculated what the mass loss rate 
would be at the current positions of o ur photoevaporat- 
ing ta ils, also based on the formula of IRichling fc Yorki 
(| 19980 . We derived values 1-5 x lO-'^Mgyr-^ which are 
a couple of orders of magnitude higher than the current 
mass loss rates (Table 4). This is not a surprise since 
most of the material is already gone from the disk and 
also the mechanism that removes the dust from the disk 
is entirely different from the classical photoevaporation 
where the flow is gas rich. 

It is interesting to compare the mass loss rates for these 
objects with that for the Vega debris disk, where lSu et al.l 
(|2005f ) find an anomalously high loss rate. Taking their 
estimated total debris mass of « 3 x lO-^M® in grains 



between 1 and 46 /zm in radius and their dwell time of « 
1000 years for these grains near the star, the mass loss is 
« IQ-i^MQ/yr. The evaporating disks are losing grain 
mass at similar or higher rates (Table 4). 

4.2. Origin of "Comet" Grains 

There are many theoretical treatments of pho- 
toevaporatio n of protoplanetary disks by hot 
stars (e.g. iJohnstone et al.l Il998t IRichling fc Yorkd 



2000t iScallv fc Clarkd l20ol llVIatsuvama et al.l l200a 



Adams et al.l 12004 iThroop fc Ballvl |2005[ ) . In general, 
these treatments agree qualitatively t hat the process 
requires lO'^ fe.g. iThroop fc Ballvl l2"005^ to lO"^ yrs (e.g. 
lAdams et al.ll20(M ) to shrink to a core d isk of size 5 — 10 
AU that is resistant to furth er erosion (jThroop fc Ballvl 
120051 IScallv fcClaTy 120011 ). Our results suggest that 
the time scales toward the short end of this range are 
likely to be applicable. The evaporation is powered by 
dust absorption of UV photons, followed by gas heating 
primarily by the photoelectric effect, with the resulting 
outflow of hot gas carrying along small, sub-micron 
sized dust grains. The larger dust gr ains will be left 
in a r elatively gas- free environment. IThroop fc Ballvl 
(|2005f ) show that these grains will be in an unstable 
configuration that will collapse to the disk mid-plain 
and that planet growth could be triggered as a result. 

There are other consequences. With the gas removed, 
its damping effects on grain motions (e.g., circularization 
of orbits) are also removed. Grains will collide and initi- 
ate coUisional cascades, that lead to a highly elevated 
rate of production of small, second-generation grains. 
These grains can be ejected by photon pressure, either 
due to photons from the central star or, more likely, those 
from the 0-star. This mechanism is a natural way to pro- 
duce the large numbers of small grains required in our 
models of the infrared output of the tails of these sources. 
An additional driver for a large grain collision rate is the 
effect of photon pressure from the 0-star perturbing the 
orbits of small grains that nevertheless remain gravita- 
tionally bound to the low mass star. 

The conditions for photon-pressure blow-out of grains 
in the disk around the low mass star can be derived as 
follows. Let (3 be the ratio of photon to gravitational 
force on a grain in the vicinity of the O star by itself. 
Photon-pressure-driven blow-out occurs for f3 > 0.5. The 
equivalent ratio of forces for a grain near the low mass 
star must include its gravity, resulting in a reduction of (3 
approximately in proportion to the ratio of gravitational 
forces from the two stars at the position of the grain. It 
is easily shown that the modified ratio has a value of 0.5 
at a distance from the low mass (mass AI2 ) star given 

by 



"/d 



/ 2/3 Ml 
V M2 
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(5) 



where d is its distance from the high mass star (mass 
Ml). The term 7 is o f order 1 and de p ends on where the 
grain is in its orbit. iLamv fc PerrinI (1199 7*) have deter- 
mined values of f3 for grains of various compositions and 
sizes around a number of stars. The values are roughly 
similar for a given grain size; grain composition enters 
as a secondary parameter. We will use typical values in 
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TABLE 4 
Physical properties of the photoevaporating sources 
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[pc] 


[yr] 
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NGC 2244 


0.22 
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0.016 
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X 


10- 


-11 _ 
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-11 


IC 1396 
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lO'-^ 
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-'' 


0.012 


5 X 


10- 


-s 


2.4 


X 


10 


-12 _ 
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10- 


-11 


NGC 2264 


0.12 


0.32 


lO'-^ 




3.5 


1.47 X 


10- 


'' 


0.013 


4 X 


10- 


-V 


3.4 


X 


10 


-11 _ 


3.4 X 


10- 


-lU 



* no reliable radial velocity dispersion information available for this cluster so we used a typical value for clusters 
of similar age 



the following discussion. Taking the values for the 09.5 
V star C Oph, grains of radius between 0.01 and 1 /im 
have P « 1000 - 10000. We take a typical distance be- 
tween our high and low mass stars to be d = 0.2 pc, and 
also assume 40 Mq for the high mass star and 0.6Mq for 
the low mass one. We then find that the 0-star photon 
pressure on such small grains becomes dominant at a dis- 
tance of w 110 AU from the low mass star. Small grains 
well inside this radius will be held in orbit by the grav- 
itational field of the low-mass star. For larger grains, (3 
drops rapidly; typically it is < 50 for 10 /im radius grains. 
The corresponding distance for them where photon pres- 
sure from the O star overcomes the gravity of the low- 
mass one is w 500 AU, that is, ou tside the boundary of 
a typical protoplanetary disk fe.g.. [Andrews fc Willianii 
[2005). Therefore, the larger grains are likely to remain 
bound to the low mass star throughout the disk, where 
they will continue to be ground down by collisions. 

That is, the grains making up the "comet tail" prob- 
ably originate from coUisional cascades in the outer re- 
gions of the circumstcUar disk. These regions have pre- 
viously been cleared of gas by photoevaporation and the 
remaining solid particles are settling toward the disk mid- 
plain and starti n g to assemble into larger bodies (see 
iThroop fc Ballyl ( 2005( 1 for details), leading to genera- 
tion of these grains in a vigorous episode of coUisional 
cascades. However, given typical sizes of protoplanetary 
disks, this process is likely to be effective only up to 
a maximum grain size of about a micron radius. The 
tails are composed primarily of grains of ~ 0.01 ^m to 
« 1 /im in radius (where the lower limit is set by the 
model fits to their surface brightness profiles; Balog et 
al. 2006) . We have modifie d the upper mass limits de- 
rived by iBalog etldl (|2006f ) to reflect the approximate 
grain size limit of 1 /im. This approximate limit results 
from the tendency of larger grains to remain gravitation- 
ally bound to the low-mass star. 

4.3. Timescales for Disk Evaporation 

From the mass loss rates and typical disk masses in 
Section 4.1, we can estimate an approximate timescale of 
10^ — 10® yrs for which this phenomenon is visible. There- 
fore, the phenomenon might be quite common, however 



we see only three cases in our GTO survey of about 20 
O stars (there are additional cases around three 0-stars 
in the W5 region, Koenig et al. in preparation) suggest- 
ing that this is probably a short lived rather rare phe- 
nomenon (only about 1/4 of the observed 0-stars have 
cometary structures in their neighborhood). 

5. CONCLUSIONS 

We present HST/NICMOS Paa images and IRS spec- 
tra of cometary structures detected in Spitzer/MIPS 24 
/im images. We estimate an upper limit to the amount of 
gas in the comets' disk and tail and find that the gas-to- 
dust mass ratio is much lower than the value observed in 
the ISM: the tails are essentially gas free. Using this new 
observation we are able to constrain the flow velocity and 
thus the mass loss rate. The new mass loss rates allow 
us to estimate the timescale on which the phenomenon 
occurs (10^ — 10^ yr). The short timescale favors photo- 
evaporation models predicting quick removal of gas from 
the outer parts of the disk. The "comet tails" are pro- 
duced from the outer regions of the disks, where larger 
grains collide at an elevated rate generating second gen- 
eration dust. These small grains are then ejected due to 
photon pressure from the nearby 0-star. The SED of 
the sources shows excess emission between 3 and 8 /im, 
in agreement with the IRS low resolution spectra. This 
emission indicates that there is an inner disk that sur- 
vives the photoevaporation process for longer than 10^ 
years, as predicted by the photoevaporation models. 
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